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Abstract 

Relativistic spin 1/2, as represented by Susskind's 1977 discretiza- 
tion of the Dirac equation on a spatial lattice, is shown to follow from 
basic, not typically relativistic but essentially quantum theoretic as- 
sumptions: that position eigenstates propagate to nearest neighbours 
while respecting lattice symmetries modulo gauge transformations. 
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1 Introduction 



The mutual consistency of quantum mechanics and special relativity has 
remained a nontrivial issue, particularly with respect to locality |^ and 
quantum measurements [0], but also in describing a free particle. Quan- 
tum mechanics, on the one-particle level, makes a fundamental distinction 
between the roles of space and time which, it seems, can be overcome only 
in the framework of quantum field theory. But even in quantum field theory, 
in the usual line of argument, relativity is something that needs to be en- 
forced. Moreover, if it is enforced by postulating unitary representations of 
the Poincare group, there appear theoretical possibilities which are never 
observed in Nature, such as continuous spinQ or tachyons. 

In this paper I suggest an ahistorical route to relativistic quantum me- 
chanics as represented by the Dirac equation. I derive the equation, and 
with it the Lorentz invariance, from seemingly "non" -relativistic quantum 
theory. In fact, Dirac himself came to the conclusion in the 1950s 0] that 
the Michelson-Morley experiment, in view of symmetries present in quantum 
but not classical mechanics, had been overinterpreted as a support of special 
relativity. As for a derivation of Lorentz invariance from a "mechanism", 
there is a well-known precedent: Maxwell in a balance of working hy- 
pothesis and actual belief [^, utilized mechanistic ideas of electromagnetic 
fields which did not enforce relativity but got it right automatically. More 
recently, in constructing cellular automata, Bialynicki-Birula noted that 
an automaton simulating the Weyl equation would require only very general 
conditions: a two-component wave function, an evolution that is linear and 
unitary, and (a vague remnant of relativity) that a wave function constant 
in space be also constant in time. 

By the technicalities used (not by the direction of argument) the present 
paper is based on a discretization of the Dirac equation devised by Susskind 
in 1977 [|]. The issue was to compensate for the doubling of the degrees of 
freedom encountered in replacing a derivative with an antihermitian differ- 
ence 

A/ _ f{x + a)-f{x-a) 
Ax 2a 

A zero difference function, for example, is not only obtained from / = const 
but also from an alternating constant (—1)^/"^ on the lattice sites. Susskind 

^Continuous spin arises if, in terms of induced representations, the little group of a 
light-like four-momentum is represented non-trivially in all of its components. 
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showed that spinorial degrees of freedom can be consistently assigned to 
different sites on a 3-dimensional lattice, thus thinning out the degeneracy 
of energy levels on a given lattice by a factor of four. The discretized Dirac 
equation resulting in this way is, in case of zero mass. 



where tp{x,y, z,t) is a one-component wave function. Thus spin 1/2, usually 
thought of as "internal" to a point particle, can be encoded in a spatial 
arrangement of hopping amplitudes for a particle without internal structure. 

In fact, equation (|l|) is the unique consequence of basic, not typically 
relativistic assumptions on the propagation of quantum particles living on 
the sites (as opposed to links or plaquettes) of a cubic spatial lattice: 

• Locality: immediate propagation to nearest neighbours only 

• Lattice symmetries are realised modulo gauge transformations 

Also, time evolution will be assumed to be linear and unitary. There will 
remain two kinematical options, one of which will be discarded because it 
is infinitely slower than the other. Remarkably, the slow option is the one 
that would realise lattice symmetries in a strict sense, without accompanying 
gauge transformations. 

In Section 2 the assumptions are specified; they include a general equa- 
tion for linear, unitary propagation as it was already proposed by this author 
0. In Section 3, the assumed invariances of the equation of propagation are 
evaluated, and equation (|1|) is derived. As for introducing particle mass, it is 
pointed out that an alternative to the standard term due to Susskind exists 
which avoids species doubling on infinite lattices. In Section 4, I present 
my Conclusions. In the Appendix, some omissions and simplifications antic- 
ipated in Sections 2 and 3 are justified. 

2 Specifying the assumptions 
2.1 Nearest-neighbour hopping 

We assume, as in 0, that a quantum particle, initially in a position eigen- 
state, will "move" by gradually (differentiably in t) forming superpositions 



itp{x,y,z,t) 



= i{ij{x + l,y,z,t)-tlj{x-l,y,z,t)) 

+ i {^{x, y+l,z,t)- ipix, y-l,z, t)) (-1) 

+ i {^{x, y,z + l,t)- ipix, y,z-l, t)) (-1) 



(1) 



x+y 
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of nearest-neighbour eigenstates. Then a general state vector, given as a 
superposition of eigenstates with coefficients ip{s,t), will evolve according to 

i'4'{s,t) ='^K{s,n) 'ip{s + n,t) (2) 

n 

where the sum runs over nearest neighbours, and where k{s, n) are complex 
hopping amplitudes whose properties are to be determined. The sum also 
includes an on-site hopping amplitude represented by n = 0. 

Since the distances between nearest neighbours are all the same, we as- 
sume that all hopping amplitudes are of the same magnitude. With a suitable 
rescaling of time we thus assume 

n)\ = 1 for all s and all n 7^ (3) 

Unitarity of time evolution and hermiticity of the Hamiltonian will be taken 
for granted. Using the standard scalar product of wave functions, 

s 

the linear operator acting on the RHS of (0) is hermitian if and only if 

k{s, —fi) = k{s — n,n) (4) 

2.2 Invariances modulo gauge transformations 

For a free particle, the equation of motion should be "the same" at all times 
and locations, as well as after a rotation. Quantum mechanically, the arbi- 
trariness of the phases of position eigenstates allows to interpret "the same" 
as "gauge equivalent" . 

In a local gauge transformation, the wave function at each space-time 
point is multiplied by a phase factor. Thus 

V'(s,^)old = g{s,t)lp{s,t)new \9is,t) \ = 1 

In terms of the new wave function, equation (0) involves the hopping ampli- 
tudes 

K{s,n)new = 9 {s + n,t) k{s, ft) old 9 {s,ty^ nj^O (5) 

K{s,0)new = k{s , 0) old - Wis, t) g{s, t)'^ (6) 
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If is a symmetry operation on the lattice (translation or rotation) the 
hopping amplitudes of equation (Q) would in general change according to 

Our assumption is that Knew is a local gauge transform of KoM^ 

K(s,n)new = 9(3 + u) n)oid 
Expressing this entirely in terms of Hoid, and dropping the index, we have 

g{s + n) K{s,n) g{s)~^ = k(S~^s, S^'^n) (7) 
Since g depends on S we eventually write g{s, S). 

2.3 Maximal gauge fixing 

Working in a particular gauge will greatly faciliate the evaluation of sym- 
metries up to gauge transformations. Following the procedure of maximal 



gauge fixing as devised in |]10l for Hamiltonian lattice gauge theories, let us 
choose to have 

y, z,l) = 1 for all x, y, z (8) 

This is accomplished, using (H), by a gauge transformation with a suitable 
behaviour in the 1 direction: 

g{x + 1, y, z, t) = g{x, y, z, t) k{x, y, z, 

The values oi g on a plane with a constant x coordinate are still free, and 
can be used to fix k,{s, 2) on that plane. Let us choose to have 

K{0,y,z,2) = lfoT al\y,z (9) 

which requires 

g{0, y+l,z,t) = g{0, y, z, t) k{0, y, z, 2)^^l 

Finally, the values of g along the line x = y = can be chosen so that 

k(0,0,z,3) = 1 (10) 

Any further gauge transformation g that is not constant throughout the 
lattice will destroy at least one of the conditions (|8D- (p^ . 
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3 Propagation on a simple cubic lattice 

3.1 Symmetries used 

For a translation by a vector a we have 

S^^n = n S^^s = s — d (11) 

For a rotation by 90° about the 1 axis, 
S-^ 1 = 1 

S-'2 = -3 S-\x,y,z) = {x,z,-y) (12) 

S-^ 3 = 2 

For a rotation by 90° about the 3 axis, 

s-^ i = -2 

S-'2 = i S-\x,y,z) = {y,-x,z) (13) 

S'^ 3 = 3 

3.2 Determining the hopping amphtudes 
3.2.1 Evaluating translations 

Let S in equation (^ be a translation as specified in (0). Thus 

g{s + n,a)K{s,n)g{s,a)^^ = K{s — d,n) n^O (14) 

Putting n = i in ( p^ and using (|]) we see that g{s, a) must be independent 
of the coordinate x, 

g{x,y,z,d) = g{y,z,d) for all a 
Now putting n = 2 and a = 1 in ([T^), we find 

/t(x - 1, z, 2) = k{x, y, z, 2) (^g{y + 1, z, 1) ^(y, 2;, i)"^) 
Solving the recursion in x and using gauge condition (H), we have 

Kix, y, 2) = e^^"(^'^) where e^"(^-^) = giy, z, 1) (7(2/ + 1, z, ly' 
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In fact, a must be independent of y and z since by re-inserting the last equa- 
tion in (|l^ (with n = 2) and considering a = 2, 3 we encounter an x depen- 
dence of exp{ixa{y, z)) on the LHS and exp{ixa{y—l, z)) or exp{ixa{y, z—1)), 
respectively, on the RHS. Hence, a{y — l,z) = a{y, z — 1) = a{y, z) so that 

/€(a;,y,^,2)=e^"^- (15) 

Using (|1^), equation ([T^) with ri = 2, a = 1, 2 can now be read as a recursion 
relation determining the gauge transformations g{s, 1) and g{s, 2) up to their 
values at X = ?/ = 0. We find 

g{y,z,i) = ^(0,^,i)e^^"^ (16) 
g{y,z,2) = g{0,z,2) (17) 

To obtain restrictions on the hopping amplitude in the 3 direction, we now 
insert ([TUp and (|T^) into (^^, using n = 3 and a = 1, 2. Thus 

g{0, z + 1,1) K{x,y,z, 3) g{0,z,i)-^ = K{x-l,y,z,3) 
g{0, z + 1,2) K{x,y,z, 3) g{0,z, 2)-^ = k(x, ?/ - 1, z, 3) 

Taking into account the gauge condition (^) the recursions are readily re- 
solved, yielding 

fi:(x,|/,2,3) = e"^(") e'^^(^) (18) 



where 



e^'^(^) = g{0,z + l,i)'' g{0,z,i) 
e*7W = g{0,z + l,2)-' g{0,z,2) 



3.2.2 Evaluating unitarity 

Applying (g) to (|), (|ll), (|1|) we obtain 



K{x,y,z,-i) = 1 

K{x,y,z,~2) = e-*-^ (19) 
k{x, y, z, -3) = e-'"'^("-') e'^^'^^"-^) 
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3.2.3 Evaluating 90° rotations about the x axis 

Let S in equation be the rotation specified by (p^. Putting ri = 1, 2, 3, 

g{x+l,y,z) K{x,y,z,l) g{x,y,zY^ = K{x,z,-y,i) 

g{x,y + l,z) K{x,y,z,2) g{x,y,z)-^ = k{x, z, -y, -?>) (20) 

g{x,y,z + 1) K{x,y,z,?>) g{x,y,z)~^ = K{x,z,-y,2) 

whence, using (I),®, (^) and (|19D, 

g{x + l,y,z) g{x,y,zy^ = 1 
g{x, y + l,z) e'°" g{x, y, zy' = e-^^^"^"') ^-i-ii^y-i) 
g{x, y,z + l) e"^(") e'^^^-^'^ g{x, y, z)-' = e'"" 

By the first of these equations, g must not depend on x. Thus, in the third 
and second equation, the only x dependence occurs in the exponentials, im- 
plying 

Thus f3 = const, and there remain two possibilities, 

gi« = = ±1 (21) 

3.2.4 Evaluating 90° rotations about the z axis 

Now let S in equation (|^) be the rotation specified by ([T3|) . Putting n = 1,2, 

g{x + l,y,z) K{x,y,z,i) g{x,y,zy^ = K{y,-x,z,-2) 
g{x,y + l,z) K{x,y,z,2) g{x,y,z)-^ = K{y,-x,z,i) 

Using (|),(|ll),(|ll) we obtain 

g{x + l,y,z) g{x,y,z)-^ = e"^"^ 
g{x,y + l,z)g{x,y,z)-^ = e"^"^ 

The solution to these recursion relations is 

gix,y,z)=e-'''^y giO,0,z) 

The xy dependent factor drops out when inserted in (|^ with n = 3, leaving 

g{0, 0,z + l) k{x, y, z, 3) g{0, 0, z)'^ = K{y, -x, z, 3) 



Using (|TS|) we arrive at 



^(0, 0,z + l) e*"^ e'y^^'^ g{0, 0, z)-^ = e'^^ e'^^^") 

Considering the x and y dependences we obtain 

e'^ = e'" (22) 

Thus, using (|2T|) , the hopping amphtudes are determined up to the choice of 
a = 0, 

K{x,y,z,n) = l (23) 

or a = TT, 

K{x,y,z,l) = K{x,y,z,—1) = 1 

K{x,y,z,2) = t,{x,y,z,-2) = (-1)^ (24) 
K{x,y,z,3) = K{x,y,z,-3) = (-1)^+?^ 

3.3 Staticity of the scalar solution 

Option (^) for the hopping amphtudes would also result from postulating 
strict invariance under the lattice symmetries, as already studied in |^ . Its 
continuum limit was found to be the nonrelativistic Schrodinger equation 
of a scalar particle. However, the time scaleQ on which the wave functions 
would evolve was found to be A^/zta^, where k denotes the nearest-neighbour 
hopping amplitude, and A is a length scale of the wave function. A^/zta^ is 
also the time scale of unitary cellular automata simulating scalar particles 



ri| . In contrast, the time scale of option (^) can be seen to be X/na from 
the initial choice of scale in (|]) and its modification in (|2^) ; again, k denotes 
the nearest-neighbour hopping amplitude. If a is very small (like the Planck 
length), we obviously have 

(X/af > A/a 

so option (ESf) tends to a static (non-kinetic and, in this sense, non-particle) 



scenario relative to (|2^) . 



is an inverse time by equation (^. The lattice spacing a emerges from Taylor 
expansions of next-neighbour terms; hence, its dimension is always cancelled by that of a 
spatial derivative. The derivatives act on wave functions in the continuum limit, so they 
are independent of a. 
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3.4 Recovering the massless Dirac equation 



This section reviews standard procedure with lattice fermions. Using (p^ 
the hopping equation @ reads 

ii>{x,y,z,t) = {^/j{x + l,y,z,t) + 'ip{x-l,y,z,t)) 

+ {tPix, y + l,z,t)+ ^{x, y-l,z, t)) (-1)" 
+ {iPix, y,z + l,t)+ ^{x, y,z-l, t)) (-1)^+^ 

Equation (|l|) is recovered by the gauge transformation ipoid = i^'^^'^^'^ncw 
For the new ip, 

iij{x,y,z,t) = i{ilj{x + l,y,z,t)-ilj{x-l,y,z,t)) 

+ i {-ipix, y + l,z,t)- tpix, y-l,z, t)) (-1)^ 
+ t {^{x, y,z + l,t)- ^{x, y,z-l, t)) {-l^+y 

Since the alternating sign factors are strongly fluctuating when viewed on a 
length scale much larger than the lattice spacing a, there can be no smooth 
solution to equation (|^). However, the equation is solved by a superposition 
of wave functions of the form 

ipooix, y,z,t) iljoi{x,y,z,t){-l)y , . 

ijio{x,y,z,t){-ir ^i,{x,y,z,t){-l)y+- ^ ^ 

where ipAB{x,y, z,t) is assumed to be smooth in the sense that it varies 
from a lattice site to the next in 0{a) at most. As suggested by the double 
index, the space of solutions is a tensor product. Multiplication by (— 1)^, 
for example, interchanges the presence/absence of that factor in the wave 
function, hence it is represented by the Pauli matrix ai acting on the first 
index, and by ui ® 1 acting on both indices. Similarly, differentiating along 
the X direction gives an extra minus sign depending on whether the factor 
(—1)^ is present or absent; this corresponds to the action of ® 1. Thus the 
right-hand side of (|Tp combines matrix factors and spatial differences into 

(X3 (E) 1 i^x + o"! eg) (Ta i^y + o"i (8> o"! i^z (26) 

The tensor products are readily seen to satisfy the algebraic relations of the 
Dirac a matrices. The difference operations asymptotically tend to 2ad/dx, 
2ad/dy, ladjdz in the continuum limit a — > 0. We may absorb the factor 
of 2a in a redefinition of the time parameter, thus recovering the massless 
Dirac equation 

dxh dih . , 
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3.5 Mass terms 



In spatial continuum, massless Dirac particles have chiral symmetry. This 
symmetry gets broken if a mass term is introduced. Therefore it is not 
unsatisfactory to find that mass terms on a lattice may require the breaking 
of a lattice symmetry. Susskind's mass term is an on-site hopping amplitude 
^(^—lY+y+^ so that equation (|I|) becomes 

y, z, t) = (lA, + i{-lT/\y + «(-l)"+^A, + z, t) 

The last term breaks the invariance (modulo gauge transformations) under 
translations by one lattice unit, while invariance under translations by two 
units is preserved. Susskind's mass term requires a doubling of the dimension 
of the space of solutions, since the functions (|25|) need to be complemented by 
analogous functions with an extra factor of (—1)^. This expands the tensor 
products (p6D of the Hamiltonian to 



0-3 (g) 1 (g) 1 iAx + CTi (8 (T3 (8 1 iAy + cji ai a-^ iA^ + fieri (g) ai Ui 



The enlarged space of states is also recovered by acting on (|25|) with a sym 



metry of equation Susskind's version of the parity operation 

^P{x,y,z,t) -y,-z,t) 

In Monte Carlo simulations, which can only use lattices with a finite number 
of sites, there is no natural distinction between smooth and strongly fluc- 
tuating wave functions. Thus the degeneracy related to the above parity 
operation can only be suppressed at the expense of some arbitrariness. 

On infinite lattices, however, an extra factor of (—1)^ to the functions 
p5| ) does make a difference. It may therefore be of interest to note that 
Susskind's term is not the only possibility of introducing mass. For example, 
we may allow for a variation of the magnitude of the hopping amplitude in 
the X direction so as to violate (H) while keeping (||), 

k{x, y,z,i) = i + iji (—1)^ K,{x, y, z, —1) = —i + iji (—1)^ 

Physically this would correspond to an alternating variation of the lattice 
spacing. The additional term in hopping equation (|l|) is 

ifxi-lf {^{x + l,y,z,t) + tp{x - l,y,z,t)) 
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Repeating the arguments that lead to (^) we obtain the operator 

(73 (g) 1 iAx + (7i ® 0-3 iAy + (Ji (g) (Ji iAz + 2i^{oia-i) ® 1 

Since icJicr-i = cr2 and since (12 ® 1 anticommutes with the first three tensor 
products we recover the Dirac matrices in the form 

0-3 (g) 1 = ai (Ji (g) 0-3 = ^2 CTi (g) (Ji = ^3 (T2 (g) 1 = /9 

4 Conclusions 

We have derived Susskind's discretization of the Dirac equation from assump- 
tions which apparently do not anticipate special relativity. While time was 
assumed to run continuously, spatial coordinates were confined to a lattice 
(reminiscent of a stack of particle detectors). Locality, too, was imposed in 
an unrelativistic sense, assuming that propagation from some position will, 
within a short interval of time, reach the nearest neighbours only. The in- 
trinsically quantum-mechanical assumption was that the amplitudes of prop- 
agation will respect the symmetries of the lattice to the extent they have to 
in quantum theory, namely up to phase shifts of position eigenstates. 

On the basis of these assumptions, a non-relativistic and a relativistic 



option appeared at the same stage in section p.3| . As it happened, this 
was simultaneously the alternative between strict symmetry and symmetry 
modulo gauge transformations. It was a matter of kinematical speed, rather 
than principle, that the non-relativistic option was discarded. 

What insight do we gain by this route to the Dirac equation? I think 
it explains the preferred role of spin 1/2 in the Standard Model, since no 
internal structure of a particle (of a kind living on lattice sites) was assumed, 
and yet the Dirac equation resulted. In particular, in the case of zero mass 
where continuous spin is a possibility consistent with Poincare invariance, it 
was just the massless version of the Dirac equation which emerged. More 
generally, I think the unity rather than mere consistency of special rela- 
tivity and quantum theory — even in a "non" -relativistic formulation of the 
latter — is emerging here. Finally, gauge transformations turn out to be as 
fundamental to the propagation of free particles as to particle interactions. 

It would be interesting to determine the "internal" degrees of freedom on 
lattices with other than simple cubic structure, especially with some of the 



infinitely many close-pack structures |12 
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Appendix 

Time-dependence ruled out 

To simplify the notation, hopping amplitudes were so far considered as func- 
tions of the spatial coordinates only. Here we show that in the gauge we were 
using there is, in fact, no other possibility consistent with time-translation 
invariance modulo gauge transformations. 
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Let g{s,t) be the gauge transformation accomplishing the shift S of the 
hopping amphtudes by a time St. Equation then reads 

g{s + n, t) k{s, n, t) g(s, t)'-^ = k{s, n,t — 5t) (28) 

Putting n = i and using our gauge condition (||), which holds at all times, 
we obtain g{x + 1, y, z, t) = g{x, y, z, t) and hence 

g{x,y,z,t) = g{y,z,t) 

Now putting n = 2 and x = in (|28|) and using gauge condition (|) we 
moreover obtain 

g{y,z,t) = g{z,t) 

Finally, putting n = 3 and x = y = and using gauge condition (|TUp we see 
that the gauge factor can only be a function of the time parameter: 

g{z,t) = g{t) 

But then g has no effect at all in equation (|28|) , and the nearest-neighbour 
hopping amplitudes must be strictly invariant under a time shift: 

k{s, n, t) = k{s, n,t — St) 



On-site hopping gauged away 

On-site hopping amplitudes consistent with the spatial symmetries (modulo 
gauge transformations) can always be gauged away, as we now show. Thus 



it was justified to omit them in the previous sections (excluding section |3l5 
where translational symmetry was partially broken). 

The analogue of equation (^, using gauge transformation (|), would be 

k{s, 0, t) - ig{s, t)g{s, t)-^ = k{S^^s, 0, t) (29) 

Let us reconsider the gauge factors which accomplished spatial translations 
of the nearest-neighbour amplitudes (section p. 2.1 ), taking into account the 



final expressions of the amplitudes as given by (^). Equation (|T^) with n put 
equal to 1, 2, 3 can then be read as a recursion relation determining the x, y, z 
dependence of g{x, y, z, a) for a given shift vector a. The arguments leading 
to those expressions did anticipate that k{x, y, z, h) would be independent of 
time, but this was justified in the previous section. 
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Those recursion relations implicit in (P^) determine g{x,y, z,a), at each 
instant of time, up to a global phase factor. That factor could depend on t, 
allowing for an expression of the form 

g{x, y, z, t, a) = go{x, y, z, a) ■ h{t, a) 



Inserting this in (p9D the time-independent go drops out, leaving 

k{s, 0, t) — ih{t, d)h(t, d)^^ = k(s — a, 0, t) 

In particular, putting a = 1, 2, 3 and defining 

c = (ci,C2,C3) c„(t) = ih{t,n)h{t,h)~'^ 

we see that on-site hopping amplitudes can at most take the form 

k{s,0, t) = k{0,0, t) + s- c{t) (30) 

But c{t) must vanish due to rotational symmetry (modulo gauge transforma- 
tions) by an argument similar to the above for translations. Reconsidering 
section p.2.3| we find that equations ( [20| ) in conjunction with (^) determine 
the gauge factor up to a time-dependent global factor, 

g{x, y, z, t, R-,) = goix, y, z, R^) ■ h{t, R^) 

In equation (pQ]), again, the time-independent g^ drops out, leaving 

k(s, 0, t) - ih{t, R^)h{t, R^)-^ = k{R-^s, 0, t) 

Taken at the origin s = the equation implies that h must vanish. Hence, 
n{s, 0, t) must be strictly invariant under the 90° rotation about the x axis, 
which implies that c(t) can at most have an x component. This latter possi- 
bility can finally be ruled out by reconsidering the 90° rotation about the z 
axis as in section ^.2.4 . 



The remaining term of (^) can be removed by a gauge transformation 
only dependent on time, satisfying g{t) = —iK(0,0,t)g(t). 
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